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ABSTRACT 

We investigate the dynamical basis of the classic empirical models (specifically, Sersic-Einasto 
and generalized NFW) that are widely used to describe the distributions of collisionless matter in 
galaxies. We submit that such a basis is provided by our a-profiles, shown to constitute solutions 
of the Jeans dynamical equilibrium with physical boundary conditions. We show how to set the 
parameters of the empirical in terms of the dynamical models; we find the empirical models, 
and specifically Sersic-Einasto, to constitute a simple and close approximation to the dynamical 
models. Finally, we discuss how these provide an useful baseline for assessing the impact of the 
small-scale dynamics that may modulate the density slope in the central galaxy regions. 

Subject headings: dark matter — galaxies: halos — galaxies: structure — methods: analytical 



1. Introduction 

The classic Sersic (1963) models met a wide 
and lasting success as empirical representations of 
the projected (2-dimensional) light distributions 
in spheroidal galaxies (for a review, see Kormendy 
et al. 2009). Einasto (1965) developed and used 
a similar shape to describe in simple terms 3- 
dimensional stellar mass profiles. 

On the other hand, recent extensive iV-body 
simulations (e.g., Navarro et al. 2004; Merritt et 
al. 2005; Gao et al. 2008; Stadel et al. 2009; 
Navarro et al. 2010) indicate that the Sersic and 
Einasto functional forms also provide good pat- 
terns to represent the spherically-averaged mass 
distributions in dark matter (DM) halos ranging 
from galaxies to galaxy clusters. These apply 
at levels comparable to, or even better than the 
popular NFW formula (Navarro, Frenk & White 
1997). 

Still, no agreed understanding is available to 
explain the value in both the real and the virtual 
world of the Sersic and Einasto representations 
(see discussions by Graham et al. 2006; Kormendy 
et al. 2009). Can we identify the underlying as- 



trophysical basis? 

2. Empirical models 

Before addressing the issue, we note that these 
models belong to two main families: generalized 
NFW (see Hernquist 1990; Zhao 1996; Widrow 
2000; hereafter gNFW) and Sersic-Einasto (see 
Graham et al. 2006; Merritt et al. 2006; Prug- 
niel & Simien 1997; hereafter SE). 

2.1. Density runs 

The density runs of the SE family may be rep- 
resented in the form 
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Here, quantities are normalized to their value at 
r_ 2 , the reference radius where the logarithmic 
slope 7 = — dlogp/dlogr takes on the value 2; 
typically, in nearby elliptical galaxies r_2 corre- 
sponds to sizes of order 10 kpc, a few times the 
half-light radius R e . 

The parameters t and r\ describe the inner slope 
and the middle curvature of the density run, re- 
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spectively. The original Einasto profile belongs to 
this family, and is obtained when r = 0. Note, 
however, that by deprojecting from the plane of 
the sky a Sersic 2-dimensional run e~ s 1 with in- 
dex n « 3 — 4 (suited for normal ellipticals, see Ko- 
rmendy et al. 2009) produces a cuspy inner run as 
in Eq. (1) with r ~ 1 - 1.19/2n + 0.22/4n 2 w 0.8 
significantly different from and less than 1, as 
shown by Prugniel & Simien (1997). 

On the DM side, recent simulations (see Gao 
et al. 2008; Stadel et al. 2009; Navarro et al. 
2010) only provide an upper bound r < 0.9 for 
the inner slope. When the original Einasto profile 
(with t = 0) is adopted, the best-fit to simulated 
DM halos obtains for rj m 0.2; we will come back 
to this value later on. 

In turn, the density runs of the gNFW family 
may be written in the forrr0: 

P{r)=r [l + ^fn) ' W = ' 2-r-^ ; 

(2) 

the parameters r, rj, and £ describe the central 
slope, the middle curvature, and the outer de- 
cline of the density run, respectively. Note that 
familiar empirical profiles are recovered for spe- 
cific values of the triple (r, t], £); e.g., Plummer's 
(1911) corresponds to (0,2,2.5), Jaffe's (1983) to 
(2,1,2), Hernquist's (1990) to (1,1,3), and NFW 
to (1,1,2). 

2.2. Toward a single family 

The main apparent difference between SE and 
gNFW is constituted by the former's exponen- 
tial decline vs. the latter's powerlaw falloff p oc 
r -{T+ni) f or large r. 

On the other hand, Eq. (2) is to be considered 
for large values of £ anyway, since a steep density 
run in the halo outskirts is indicated by observa- 
tions of light distribution in spheroidal galaxies 
(other than cDs, see Kormendy et al. 2009), and 
of DM distributions from weak lensing in galaxies 
and galaxy clusters (e.g., Broadhurst et al. 2008; 
Oguri et al. 2009; Newman et al. 2009). 

The circumstance is easily translated into the 
formal statement that the gNFW family converges 



In the literature these runs are sometimes referred to as 
a/37-models, and equivalently defined via the parameters 
7 = t, a = I/77, = t + v£- 



to the SE for large £. This is seen on recasting pf T 
from Eq. (2) in exponential form, to read 

e £ln[(l+w)/(l+wr>)] ^ e iw(l~r") ^ e ~u(r"-l) 

(3) 

for approximating the middle and last terms we 
have used the circumstance that £ 3> 1 implies 
w >• 1 and so £ w ~ (2 — r)/r) = u applies. Thus 
the two families in Eqs. (1) and (2) actually be- 
come one in this limit. 

Thus in the following we focus mainly on the 
SE family, and proceed to discuss its dynamical 
basis in terms of the Jeans equation. 

3. The Dynamical model 

The dynamical model of DM halos hinges upon 
the radial Jeans equation that expresses the self- 
gravitating, equilibrium of collisionless matter (see 
Binney & Tremaine 2008). The Jeans equation 
reads 

I djpo*) _ GM{< r) 20 0* 
p ar r z r 

in terms of the density p(r), the related cumulative 
mass M(< r) = 4-7T fldx x 2 p(x), and the radial 
velocity dispersion of(r). The last term on the 
r.h.s. describes the effects of anisotropic random 
velocities via the standard Binney (1978) param- 
eter P = 1 - <Tg/(T 2 . 

Note that the Jeans equation is designed to de- 
scribe a (quasi-)static equilibrium, away from ex- 
treme major merger events like is the case with 
the Bullet Cluster (see Clowe et al. 2006). But 
even in relaxed conditions, solving Jeans requires 
an 'equation of state', i.e., a functional relation ex- 
pressing the DM pressure p of in terms of density 
(and possibly radius) only. 

3.1. Equation of state 

In seeking for such a relation, one can make 
contact with the classic theory of the non-linear 
collapse for DM halos in an expanding Universe; 
here self-similar arguments play the role of a 
pivotal pattern (see Fillmore & Goldreich 1983; 
Bertschinger 1985; Taylor & Navarro 2001). This 
still applies to modern views of the halo develop- 
ment (e.g., Mo & Mao 2004; Lu et al. 2006; Li et 
al. 2007; Lapi & Cavaliere 2009a,b; Fakhouri et al. 
2010), that comprise two stages: an early collapse 
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builds up the halo main body via a few major 
merger events and sets its phase-space structure 
by dynamical relaxation of DM particle orbits; 
this tails off into a secular development of the 
outskirts by smooth accretion and minor mergers. 

The essence of the macroscopic equilibrium is 
conveyed by the self-similar scaling a 2 oc M/r 
adding to the geometric relation p oc M/r 3 . The 
macroscopic import of the halo phase-space struc- 
ture is conveyed by combining these two quanti- 
ties into the 'phase-space density' p/af, or equiv- 
alently into the functional K(r) = a 2 /p 2 / 3 often 
referred to as DM 'entropy' (sec Bcrtschinger 1985; 
Taylor & Navarro 2001). For the latter quantity, 
one easily derives the scaling K(r) oc rM 1 ' 3 im- 
plying 

K{r) oc r a ; (5) 

whence one expects a slope a slightly exceeding 
unity. 

To focus the values of a, in Lapi & Cavaliere 
(2009a, b) we have developed a full semianalytic 
treatment of the halo growth in the standard ac- 
celerating Universe (see Komatsu et al. 2010). We 
found constant values of a, that fall within the 
narrow range 1.25 — 1.3; on average, such values 
grow weakly with the mass of the halo body, from 
galaxies to rich clusters. 

The halo development process has been probed, 
and the two-stage view confirmed by intensive TV- 
body simulations (e.g., Zhao et al. 2003; Wechsler 
et al. 2006; Diemand et al. 2007; Hoffman et al. 
2007; Schmidt et al. 2008; Ascasibar & Gottlober 
2008; Vass et al. 2009; Genel et al. 2010; Wang et 
al. 2010). These also confirm that: (i) a (quasi- 
static macroscopic equilibrium is attained at the 
end of the fast collapse, and is retained during the 
subsequent stage of secular, smooth mass addi- 
tion; (ii) a persistent feature of such an equilibrium 
is constituted by powerlaw correlations holding in 
the form o 2 jp 2 l 3 oc r a , although it is still widely 
debated whether the radial or the total velocity 
dispersion best applies (see also the discussions by 
Schmidt et al. 2008 and by Navarro et al. 2010). 

In building up our dynamical models we fo- 
cus on the quantity K = a 2 /p 2 ^ 3 oc r a that in- 
volves the radial dispersion a 2 (see also Dehnen 
& McLaughlin 2005). Operationally, this provides 
a direct expression for the radial pressure term 
pa 2 r = K p 5 / 3 oc r a p 5 / 3 in the Jeans Eq. (4); 



anisotropics are accounted for by the last term on 
the r.h.s., as discussed in § 3.3 below. 

3.2. The DM a-profiles 

In terms of K(r) oc r a , the Jeans equation may 
be recast into the compact form 

with 7 = — dlogp/dlogr representing the loga- 
rithmic density slope and v 2 = GM(< r)/r the 
circular velocity. Remarkably, by double differen- 
tiation this integro-differential equation for p(r) 
reduces to a handy 2 nd order differential equation 
for 7 (see Austin et al. 2005; Dehnen & McLaugh- 
lin 2005). 

Tackling first the isotropic case j3 — 0, we re- 
call that the solution space of Eq. (6) spans the 
range a < 35/27 = 1.296: the specific solution 
for the upper bound, and the behaviors of the 
others have been analytically investigated by Tay- 
lor & Navarro (2001), Austin et al. (2005) and 
Dehnen & McLaughlin (2005). In Lapi & Cava- 
liere (2009a) , we explicitly derived the solutions in 
the full range a = 1.25 — 1.296, that are marked by 
a monotonically decreasing run and satisfy physi- 
cal boundary conditions: a finite central pressure 
or energy density (equivalent to a round minimum 
of the gravitational potential); a steep outer run 
implying a finite and rapidly converging (hence a 
definite) overall mass. We dubbed a-profiles these 
physical solutions. 

We shall use the following basic features of the 
latter. In the halo body at the point ro the a- 
profile is tangent to the pure powerlaw solution 
p oc r~ 70 of the Jeans equation; there v 2 /<r 2 oc 
r 2-7o/3-a a pp ueS) to imply from Eq. (6) 

70 = 6 - 3 a . (7) 

This is consistent for a = 1.25 with the self-similar 
slope; as such, it qualifies to provide a univer- 
sal middle-range slope. Note that the point r 
lies in the neighborhood of the radius r_ 2 (see 
§ 2.1), specifically ro « 1.74 — 1.51 r_2 holds for 
aw 1.25-1.3. 

On the other hand, a monotonic density run 
implies the term v^/a 2 . oc r 2 t/ 3 ~ a to vanish at 
the center; this results in an inner powerlaw p oc 
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r 7a with the slope 



7a 



3 

5^ 



(8) 



This differs from zero as long as the entropy run 
grows from the center following K oc r a with a > 
0. 

Finally, a finite mass implies v 2 /a 2 oc r - 1 + 2 7/3-c 
to hold in the outskirts, so as to yield a typical 
outer decline p oc r~ 7b with slope 



76 = - (1 + a) . 



(9) 



This exceeds the value 3, and so constitutes the 
hallmark of a rapidly saturating mass; the circum- 
stance makes less compelling here the role of a 
virial boundary. 

Thus, compared to NFW the inner slope of 
the dynamical model is considerably flatter and 
the outer slope steeper; compared to the original 
Einasto profile, the main difference occurs in the 
inner regions where the dynamical model is (mod- 
erately) steeper. 

3.3. Anisotropy 

It is clear from Eq. (6) that anisotropics will 
steepen the density run for positive f3, and flatten 
it for negative (3. The latter condition is expected 
to prevail in the inner region, where tangential 
components develop from the angular momentum 
barrier (Nusser 2001; Lu et al. 2006). Moving 
outwards, radial motions are expected to prevail, 
so raising f3 up to values around 0.5 at r « ro; 
outwards of this, j3 is expected to saturate or even 
decrease, as one enters a region increasingly popu- 
lated by DM particles on eccentric orbits with van- 
ishing radial dispersions at their apocenters (see 
Bertschinger 1985). 

This view is supported by numerical simula- 
tions (see Austin et al. 2005; Dehnen & McLaugh- 
lin 2005; Hansen & Moore 2006; Navarro et 
al. 2010), which in detail suggest the average 
anisotropy-density slope relation 



/3(r)«/3(0)+/?' [7M-7a] 



(10) 



to hold with parameters /3(0) « -0.1, /?' « 0.2, 
and the constraint (3(r) < 0.5. Note that at all 
radii the inequality 7(7") > 2/3(r) is satisfied; this 
has been conjectured to constitute a necessary 



Table 1: Parameters of empirical models in the 
isotropic case 





1.25 


1.27 


1.29 


Einasto model (Eq. 1) 

V 


0.211 


0.182 


0.152 


SE model (Eq. 1) 

T 

V 


0.750 
0.327 


0.762 
0.287 


0.774 
0.244 


gNFW model (Eq. 2) 

T 

V 


0.750 
0.687 
3.821 


0.762 
0.579 
4.564 


0.774 
0.473 
5.624 



condition for a self-consistent spherical model with 
positive distribution function (see Ciotti & Mor- 
ganti 2010). 

In Lapi & Cavaliere (2009b) we extended the 
dynamical model to such anisotropic conditions in 
the full range a w 1.25 — 1.3, inspired by the anal- 
ysis of Dehnen & McLaughlin (2005) for the upper 
bound of a. We note that the latter is now slightly 
modified to 35/27-4/3(0)/27 w 1.31; likewise, the 
point r where 7 = 70 = 6 — 3a applies moves 
slightly inwards, so that r w 1.58 — 1.38 r_ 2 now 
holds. 

The main outcome, however, is that the density 
profile is somewhat flattened at the center relative 
to the isotropic case; the inner slope now reads 



7a = I a + jj (3(0) . 



(11) 



In particular, even a limited central anisotropy 
(corresponding to values /3(0) w —0.1) causes an 
appreciable flattening down to j a w 0.63 — 0.66 
for a w 1.25 - 1.3. 

On the other hand, we stress that such small 
phcnomenological anisotropies near the center im- 
ply the radial and the total dispersions a 2 = 
a 2 r [l-2 /3/3] to be very close. 

4. From Dynamical to Empirical Models 

Here we discuss how the parameters of the em- 
pirical profiles (see § 2) can be set based on our 
dynamical model (see § 3); in such conditions, it 
will turn out that such profiles constitute close 
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Isotropic 



Anisotropic 




Fig. 1. — Density profiles in isotropic (left panels) and anisotropic conditions (right panels) for the dynamical 
model and its approximation in terms of empirical models with the parameters derived in § 4 and summarized 
in Table 1 and 2. The lower panels highlight the corresponding logarithmic density slopes. 



Table 2: Parameters of empirical models in the 
anisotropic case 



a 


1.25 


1.27 


1.29 


Einasto model (Eq. 1) 








V 


0.259 


0.226 


0.194 


SE model (Eq. 1) 








T 


0.630 


0.642 


0.654 


V 


0.368 


0.326 


0.285 


gNFW model (Eq. 2) 








r 


0.630 


0.642 


0.654 


V 


0.808 


0.688 


0.578 




3.396 


4.018 


4.812 



approximations to the model over a wide radial 
range. 

4.1. Parameters from dynamics 

First we consider the original Einasto profile 
(r = in Eq. 1), since this has been widely used 



in the context of DM halo simulations. Here, r is 
fixed to 0, and the only free parameter is the cur- 
vature rj. This we set by requiring the logarithmic 
density slope 

j(f) = 2f v (12) 

to equal 70 at the point f ■ So we find the expres- 
sion 

, = (i3) 

logr 

that takes on values r) ~ 0.15 — 0.2, see Table 1 
and 2; remarkably, these turn out to agree with 
those derived from fits of state-of-the-art TV-body 
simulations in terms of the same Einasto density 
run, as performed by Navarro et al. (2010). 

On the other hand, the flat central slope of 
the Einasto profile is at variance with the value 
7 a = 3 a/5 given by our dynamical models based 
on Jeans; to wit, consistency between pure Einasto 
and Jeans would require at the center a flat en- 
tropy distribution and a vanishing pressure. 

Actually, the simulations quoted in § 3.1 within 
their finite mass resolution provide only an upper 
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Isotropic 



Anisotropic 




Fig. 2. — Same as Fig. 1 for the profiles of circular velocity, and for the corresponding logarithmic slopes. 



limit t < 0.9 to the central slope. This grants 
scope to the full SE family of Eq. (1). 

The latter features two parameters, the inner 
slope t and the middle curvature r\. These we set 
by requiring the logarithmic density slope 

j(f) = r + (2 — t) f (14) 

to equal 7 a for f — )• 0, and 70 at fo; so we find 

T = 7a 

(15) 

_ log[(7Q-7a)/(2-7a)] 
logfo 

Thus we predict the central slope to take on values 
r = 0.6 — 0.8 and the corresponding curvature 
parameter to take on values r\ = 0.2—0.3, see Table 
1 and 2. It will be worth fitting the outcomes of 
TV-body simulations based on these extended SE 
profiles with r > 0. 

Finally, we report the corresponding results for 
the empirical gNFW family. This features three 
parameters: inner slope r, middle curvature 77, and 
strength of the outer decline £; these we set by 



requiring the logarithmic density slope 

71 > +% [(2 - t)/(2 -r -1 1 ° j 

to equal "f a for f — > 0, 70 at fo, and 7^ for r — > 00. 
So we find 

T = 7a 

log [(70 - 7a) (2 - 7fc)/(7o - 7b) ( 2 ~ 7a)] 

V = ; ; 

log r 



rj 

The parameters so determined are listed in Ta- 
ble 1 and 2 for both the isotropic and anisotropic 
conditions. 

4.2. Results and comments 

With the parameters focused as discussed in the 
previous subsection, Fig. 1 illustrates how the em- 
pirical compare with our dynamical models. We 
plot the density run of the latter (specifically for 
the a-profile with a s» 1.25 suitable to galactic 
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halos), compared to those of the Einasto, SE, and 
gNFW models. The left and right panels refer to 
isotropic and anisotropic conditions, respectively; 
the popular NFW profile is also shown for refer- 
ence. To make comparisons easier, we plot in the 
lower panels the corresponding logarithmic density 
slopes. 

It turns out that the closest approximation to 
the dynamical model is provided by SE, which 
shares with it not only the central slope by con- 
struction, but also the body and the outer behav- 
iors. The original Einasto profile provides an ac- 
ceptable approximation in the middle and outer 
ranges, but not at the center, because of its flat- 
ness. On the other hand, the gNFW family pro- 
vides an acceptable approximation in the inner 
and middle ranges, but not in the outskirts where 
its slope is too flat. Finally, the NFW profile 
provides an acceptable approximation only in the 
middle range. 

Similar conclusions concern the profiles of cir- 
cular velocities v 2 (r) = GM(< r)/r, that are an- 
alytically dealt with in the Appendix and illus- 
trated in Fig. 2. 

We stress that the handy SE representation is 
convenient in analyzing data in several contexts, 
including: the DM particle annihilation signal ex- 
pected from the Galactic Center (see Lapi et al. 
2010a); rotation curves of dwarf and normal spi- 
ral galaxies (see Salucci et al. 2007); individual 
and statistical properties of elliptical and spiral 
galaxies (see Cook et al. 2009); strong and weak 
gravitational lensing (see Lapi et al. 2009b), cur- 
rently observed in clusters (e.g., Zitrin et al. 2010) 
and soon in massive elliptical galaxies (see discus- 
sion by Bradac et al. 2009); X-ray emission from 
the intracluster plasma (see Cavaliere et al. 2009; 
Fusco-Femiano et al. 2009, Lapi et al. 2010b). 

5. Discussion 

We first stress that the dynamical model (as 
well as its approximations in terms of empirical 
models) is in keeping with the basic features of 
standard DM, i.e., its cold and collisionless nature. 
In fact, it implies a 2 (r) — >■ for large r ^> r_ 2 , a 
behavior expected in the outskirts for cold matter 
dominating the potential well. 

At the inner end, with decreasing r we expect 
a 2 (r) to increase toward a maximum, correspond- 



ing to effective conversion of inflow kinetic into 
random energy. In fact, toward the center Jeans 
requires d log a 2 /d log r = 7 — GM(< r)/r 2 — > "f a 
to hold as the gravitational force vanishes there, 
to the effect that a 2 (r) oc r' ta -> 0. 

Concerning the collisionless nature of the DM, 
the boundary conditions at the center imply a fi- 
nite, non zero pressure (and energy density) , while 
a long collisional mean free path allows the pres- 
sure gradient dp/dr to diverge. Conversely, with a 
short mean free path A the pressure gradient can- 
not diverge on scales r > A, where a finite a 2 and 
a flatter 7 apply. In fact, weakly collisional con- 
ditions have been proposed to explain the cored 
light profiles observed in many spheroidal galaxies 
(see Ostriker 2000). 

On approaching the center of a galactic halo, 
one expects the basic dynamical model from large- 
scale Jeans equilibrium to be altered to an in- 
creasing degree by small-scale dynamics and/or 
energetics related to baryons. These processes are 
specifically related to following issues: transfer of 
energy/angular momentum from baryons to DM 
during galaxy formation; scouring baryons by the 
energy feedback from central active galactic nu- 
clei; any 'adiabatic' contraction of the baryons. 
Such issues will be briefly discussed in turn, with 
a warning that they enter increasingly debated 
grounds. 

5.1. Energy/ Angular Momentum Trans- 
fers 

Flattening of the inner density profile may be 
caused by transfer of energy and/or angular mo- 
mentum from the baryons to the DM during the 
galaxy formation process (see El-Zant et al. 2001; 
Tonini et al. 2006, Romano-Diaz et al. 2008). 

In detail, upon transfer of tangential random 
motions from the baryons to an initially isotropic 
DM structure, the density in the inner region is 
expected to behave as (Tonini et al. 2006) 

pocr -2[ 7 a+2(2- 7a )/3]/[2+(2- 7a )/3] _ ^g) 

Thus for /3 < the profile is flattened relative to 
the original 7o , down to the point of developing a 
core for /3 < - 7o /2 (2 - 7o ) « -0.3. 

However, a reliable assessment of the amount of 
angular momentum transferred from the baryons 
to the DM is still wanting, and would require 
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aimed numerical simulations of better resolution 
than presently achieved. 

5.2. Other processes on inner scales 

Less agreed processes may affect galactic scales 
r < 10 2 pc. For example, at the formation of 
a spheroid, central starbursts and a supermas- 
sive black hole may easily discharge enough energy 
(~ 10 62 erg for a black hole mass M. ~ 1O 9 M ) 
with sufficient coupling (> 1%) to blow most of 
the gaseous baryonic mass m oc r 3 ~ 7a /(3 — 7 a ) 
out of the gravitational potential well. This will 
cause an expansion of the DM and of the stellar 
distributions (see Fan ct al. 2008), that flattens 
the central slope. 

In addition, binary black hole dynamics fol- 
lowing a substantial merger may eject on longer 
timescales formed stars from radii r « 10 (M./ 
10 8 Mq) 1 /^ 3 i") pc containing an overall mass of 
a few times the black hole's, and so may cause a 
light deficit in some galaxy cores (see Merritt 2004; 
Lauer et al. 2007; Kormendy et al. 2009). A full 
discussion of the issue concerning cored vs. cusped 
ellipticals is beyond the scope of the present paper. 

5.3. Adiabatic Contraction? 

On the other hand, some steepening of the inner 
density profile may be induced by any 'adiabatic' 
contraction of the diffuse star-forming baryons 
into a disc-like structure, as proposed by Blumcn- 
thal ct al. (1986) and Mo ct al. (1998) but cur- 
rently under scrutiny, see Abadi et al. (2010). 

On the basis of the standard treatments, it is 
easily shown that in the inner region an initial 
powerlaw p{r{) oc r^ la is modified into 

p oc r- 3 /^") ; (19) 

this yields typical slopes around 0.9, steeper than 
the original -f a < 0.78 but still significantly flatter 
than 1. 

However, recent numerical simulations (see 
discussion by Abadi et al. 2010) suggest that 
the treatment of adiabatic contraction leading to 
Eq. (19) is likely to be extreme; actually, in the 
inner region the contraction is ineffective and the 
density slope hardly modified. Again, highly re- 
solved TV-body experiments are needed to clarify 
the issue. 



6. Conclusions 

We have discussed the dynamical basis of the 
Scrsic-Einasto empirical models, in terms of well- 
behaved solutions of the Jeans equation with phys- 
ical boundary conditions comprising: a finite cen- 
tral energy density, a closely self-similar body, a 
finite (definite) overall mass. 

We find the SE profile to be particularly suit- 
able to represent the general run of the dynamical 
solution. Specifically, we have discussed how to 
tune the parameters of SE in terms of the dynam- 
ical model; in such conditions, we find the former 
to constitute a simple and close approximation to 
the latter. 

The resulting SE profile shares with the dynam- 
ical model the following features: an outer steep 
decline, hence a definite overall mass; a closely self- 
similar body with slope 70 = 6 — 3a; an inner slope 
around j a = 3a/5, hence flatter than —1. The lat- 
ter slope provides an useful baseline for discussing 
alterations of the inner behavior caused by addi- 
tional baryonic processes. 

In conclusion, we submit that the dynamical 
models discussed here, namely the a-profilcs, pro- 
vide the astro-physical basis for understanding the 
empirical success of the SE profiles in fitting the 
real and the virtual observables, from galaxies to 
galaxy clusters. 
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Empirical models: circular velocities 



Here we provide analytic formulae for the circular velocity profiles related to the empirical models pre- 
sented in § 1 of the main text. The circular velocity v\ = GM(< r)/r constitutes a quantity helpful not 
only to evaluate the gravitational potential but also in the specific context of galactic rotation curves (e.g., 
Salucci et al. 2007, and references therein). 

For the SE family this comes to 
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f (3-r) fri ] 






-q ■ 




r 




'(3-r) 













(1) 



in terms of the (lower) incomplete gamma function T[a; x] = dt t a ~ x e - '. On recalling that T[a; x] ~ x a /a 
holds for x <C 1 and r[a;x] ~ T[a] holds for x » 1, one finds the asymptotic behaviors v^(r) cx f 2_r toward 
the center and v^(r) cx toward the outskirts; hence a maximum occurs at a radius r ~ r_2- 
On the other hand, for the gNFW family the circular velocity run comes to 



2-Fi 

T 


[(3-r)^ 1+ (3-r) f J 




f (3-r) f j . (3-r). " 





(2) 



in terms of hypergeometric functions 2Fi[a,b,c;x\. From their asymptotic (details may be found in 
Abramowitz & Stegun 1972) one again finds behaviors similar to the SE family. 
The full expressions are used to compute the profiles illustrated in Fig. 2. 
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